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In this work we provide a method for generating quantum circuits preparing maximally multi- 
partite entangled states using genetic programming. The presented method is faster that known 
realisations thanks to the applied fitness function and several modifications to the genetic program- 
ming schema. Moreover, we enrich the described method by the unique possibility to define an 
arbitrary structure of a system. We use the developed method to find new quantum circuits, which 
are simpler from known results. We also analyse the efficiency of generating entanglement in the 
spin chain system and in the system of complete connections. 

I. INTRODUCTION 

One of the key resources in the quantum information theory is entanglement [HE]- It is an inherent part to 
many branches of research in the field of quantum information computing. One of the most perspective applications 
of quantum computing, the quantum key distribution, utilizes the entanglement of parts of the key possessed by 
a sender and a receiver. Similarly, the experimental implementation of quantum teleportation requires sharing of 
maximally entangled states between the communicating parties. Entanglement is also necessary in applications such 
as dense coding and quantum direct communication protocols. It is considered that entanglement plays a crucial role 
in the exponential speedup of quantum algorithms [3]. 

The above applications are based on harnessing bipartite maximally entangled states. This suggests searching for 
new quantum algorithms that would make use of maximally multipartite entangled states. However, searching for 
maximally entangled states itself is a known problem El [H] . 

The main goal of this work is to provide an efficient method of searching for circuits preparing maximally multipartite 
entangled states. The algorithm described in this work is based on genetic programming (GP). We take advantage 
of efficient fitness function and introduce the modifications to the genetic programming engine that result in shorter 
runtime of the process. We also introduce the method of representing quantum circuits allowing for the reduction of 
the search space. Additionally we use narrowed sets of universal quantum gates in order to increase computational 
efficiency. Apart from improving performance of the algorithm we also take the structure of a system into consideration. 

This paper is organized as follows. In section [IT] we describe the proposed method, focusing mainly on the circuits 
representation and the fitness function. Then, in section [TTTJ we present the obtained results. 

II. DESCRIPTION OF THE METHOD 

Genetic programming is the numerical method based on evolutionary mechanisms. We decided to use such a method 
because of the two main reasons. First of all it enables to perform numerical search in complicated, mathematically 
untraceable space. On the other hand, we assume that entanglement in a quantum circuits increases in particular seg- 
ments quite independently. Genetic programming enables to exchange segments between different circuits generating 
high entanglement by the use of evolution mechanisms. 

A. General GP algorithm 

Genetic programming belongs to the family of search heuristics inspired by the mechanism of natural evolution (in 
genetic algorithms each element of a search space being candidate for a solution is encoded as a representative of a 
population). Every member of a population has its unique genetic code, which is its representation in optimization 
algorithm. In most of the cases the genetic code is a sequence of values from a finite set E of possible values of all the 
features that characterize a potential solution x 6 S™ in the search space. Searching for the optimal solution is done 
by the modification of genetic code due according to the rules of the evolution such as mutation, selection, crossover 
and inheritance. 

Mutators are the functions that change single elements of a genetic code of a population member randomly. A 
basic example of a mutator is a function that randomly changes values of a representative x at all positions with some 
non-zero probability 



M(x)i = { X A %, Probability P 

v ' \ rand(S), probability 1 - p v y 

Crossovers implement the mechanism of inheritance. This function divides parental genetic codes and create a 
new genetic code. Commonly two new codes are created at the same time from two parental codes. An example of 
such crossover is so called two point cut, where both parental codes (xi,yi) are cut into three regions and the middle 
segments are interchanged 
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where ci < c 2 are randomly chosen indices. In every iteration of the algorithm all members of the population are 
evaluated using fitness function / : S™ — > -ft which enables ordering elements. Then, using the selector function, the 
set of the best members is obtained and used to create a new generation of a population using mutation and crossover 
functions. There is a number of strategies for defining the selector function - from completely random choices to the 
deterministic choice of best representatives. 

In this work we use a solution which places somewhere in between. During every selection we begin with establishing 
a random set Rs C S of elements from population S: 

R s = rand{U C S : \U\ = n}, 

where n is a fixed number of elements in every set. Then we apply the fitness function to choose a maximally well 
fitted one from this narrowed set (a random element from a set of maximally well fitted ones) 



Selector(S) = rand{x G Rs ■ f(x) = max f(x)}, (3) 

where fix) is the fitness function. 

The strategy based on evolution mechanism makes genetic programming especially usable when parts of genetic 
code represent the features of elements of search space and can be interchanged between elements independently. In 
such case GA is expected to find the features that occur in well fitted representatives and mix them in order to find 
the best possible combination. Pseudo code representing this approach is presented in Listing 1. 

population = RandomPopulationO 
for( generationsNumber ){ 
newPopulation = [] 

forCi = 0; i<population.size()/2; i++){ 
mom = Selector (population) 
dad = Selector (population) 
(sister, brother) = Crossover (mom, dad) 
Mutator (sister) 
Mut at or (brother ) 
newPopulation. append(sister) 
newPopulation. append (brother) 

} 

population = newPopulation 

} 

Listing 1: Pseudo code representing the algorithm of genetic programming. Functions Selector, Mutator and 

Crossover work as defined in Section ITTl 

While the customization of population representation and fitness function unavoidably relies on the optimization 
problem, other parameters of genetic programming such as crossover and mutation methods are universal. When 
treating quantum circuits as the strings of integers representing quantum gates one can use various, already developed 
methods [15] . 



B. Circuits representation - population 



In order to apply the Genetic Algorithm it is necessary to define a population, a fitness function, the methods of 
crossing-over, mutation and selection. In the case of the optimization of quantum circuits generating entanglement the 
population consists of quantum circuits. The most convenient way to represent a computation in a quantum circuit 
is a sequence of quantum gates. Every quantum gate can be approximated using the gates from a set of universal 
quantum gates. In this work we use a set containing the Hadamard gate H, the R(ir/A) (called 7r/8) gate and the 
controlled-NOT gate CNOT [23] 
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However, initial experiments suggest that the i?(7r/4) gate is not necessary for generating maximally entangled 
states in small systems (i.e. with less than 7 qubits). Thus, in order to simplify the resultant state, we consider a set 
containing only H and CNOT gates when the obtained entanglement is maximal. 

One should note that using a general representation of circuits, we treat circuits that are equivalent with respect 
to permutations of qubits as completely different. In out scheme to decrease the number of equivalent circuits in the 
search space we reduce the number of available gates on a few first positions to the ones operating on first 2p qubits, 
where p is the number of a gate. Such simplification does not influence the space of feasible solutions, because it is 
not possible to act on more than 2p qubits using p gates which all act on maximum 2 qubits. 

C. Multipartite entanglement potential - fitness function 

As our goal is to find a circuit which prepares maximally multipartite entangled state, our fitness function must 
determine how entangled the resultant state is. In this case computing fitness score of a quantum circuit requires 
obtaining the output state and then the estimation of its entanglement. One should note that we do not necessarily 
need an entanglement measure to perform a successful numerical optimization. We only need a function for some 
kind of estimation of multipartite entanglement that reaches its maximum only for maximally entangled states, which 
can be characterized by a number of criteria. 

Definition 1 State p € D(X <E> y), dim(A') > dim(J^) = n is maximally entangled if and only if the reduced state is 
maximally mixed: Tr^(p) = t/n i.e. S(Ttx(p)) = n. 

As we do not have a good measure of multipartite entanglement that we can use in our computation, we try to 
find a function that estimates the entanglement quantitatively. In this work we rely on the definition of maximally 
multipartite entangled states introduced by Facchi et al. [8J. 

Definition 2 State £ C 2 is called maximally multipartite entangled state if it is maximally entangled according 
to every bipartition. 

That definition induces a construction for multipartite entanglement potential, which is based on the algebraic sum 
of the purity of all possible bipartitions of a system. Such an approach guarantees that the maximally entangled 
states are classified correctly. However, the purity function does not distinguish some states with different values of 
entanglement measures such as negativity or von Neumann entropy. For this reason in this work we use analogous 
algebraic sum of bipartite entanglement measure defined as 

41(1^)^1) = E E VN (Tr x (\i,)(m, (5) 

{x,y) 

where (X, y) is a bipartition of a system, and dim(A') > dim(y). For pure states it holds that Tr^dV'XV'l) = AA\ 
where A is the coefficients matrix of a state (i.e. \ip) = vec(A)). Thus we have: 

E VN (\^)(yj\) =E d M)M) = S(Tr y (|V)(VI)) = S(AA*) = -]>> 4 log 2 A 4 , (6) 
where S is the von Neumann entropy, and Xi are of the matrix A. 
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40 



TABLE I: The minimal number of all gates necessary to generate maximum multipartite entanglement in spin chain 

and system of complete connections 
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TABLE II: The minimal number of CNOT gates necessary to generate maximum multipartite entanglement in spin 

chain and system of complete connections 

Although we considered a number of functions, only proper bipartite measures of entanglement have met our 
restrictions. In order to select one for further computation we performed computational speed comparison (see Figure 
[TJ. Finally we decided to use algebraic sum of a measure using von Neumann Entropy over all possible bipartitions 
of a system. 

III. GENERATION OF MULTIPARTITE ENTANGLEMENT 

As an example of an application of the described method we find the minimal circuits preparing the maximally 
multipartite entangled states. We perform search both for spin chain and complete connection systems. Finally we 
compare the obtained results. The method developed in this paper is probabilistic. This makes the optimization of 
additional parameters a tedious task. In order to obtain a minimal circuit we execute an algorithm with increasing 
number of gates. We stop when a maximally entangled state is prepared in a circuit of a fixed size. When the 
hypothetical maximum of the fitness function is unreachable, we search for the maximal value of the entanglement 
potential and then optimize the size of a circuit. In this work we restrict ourselves to 8 qubits. This is caused by the 
long runtimes of the algorithm. 

The optimization is preformed simultaneously for the spin chain and the system of complete connections. By a 
completely connected system we mean a system in which we are able to act on every pair of qubits. It is equivalent 
to defining the set of all available CNOT gates as {CNOT{i]j) : i ^ j}, where CNOT(iJ) is a controlled NOT 
gate acting on i-th and j-th qubit Q . The sets of Hadamard gates and phase gates are independent of changes of 
possible qubit connections as this gate acts on one qubit only. By the spin chain system we mean a system where 
we are able to act on pairs of the nearest neighbors in a chain. In such a case the set of available CNOT gates is 
{CNOT(i; j) : \i — j\ = 1}. In each complete system we additionally analyze the minimum topology of necessary 
connections by investigating the connections graph where the vertices represent qubits of a system, and the edges 
represent connections, that is CNOT gates of a circuit connecting two qubits. 



A. Resultant states and circuits 

The upper bound of the potential of entanglement can be given explicitly if the number of all possible bipartitions is 
known. However, it is not always possible to find a pure state with the maximum value of the entanglement potential. 

We have managed to reach the hypothetical maximum in the case of 3, 5 and 6 qubits. These are the most important 
cases, because there are no concerns about the correctness of their classification by the fitness function. 3- and 5-qubit 
states obtained by us are analogous to the states presented in [20]. In the 6 qubit case the best known state was 
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TABLE III: The maximum value of the entanglement potential obtained for d — 1, .... 8 qubits and hypothetical 

upper bound 
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FIG. 1: Average computational time (log scale) required to calculate negativity and von Neumann entropy. For each 
dimension d = 1, . . . , 8 the results were obtained for a sample of 10 4 random states. 

presented in work [20]. However, that state consists of 32 non-zero coefficients. We have managed to find a state with 
16 non-zero coefficients 

|Ve) =*((|0000) - |HIl»|^ + > + (|0011> + |1100))|^->+ 
(|0101> + |1010)) |0+) + (|0110) - |lOOl»|0-», 

where |V' ± ) = |00) ± |11) and 1^) = |01) ± 1 10) . 

In the case of 4 and 7 qubits it is not known whether there is any state with hypothetical maximum of the potential 
of entanglement and in 8-qubit case it is definitely not possible. Thus we can only search for circuits that generate 
the amount of entanglement that we consider maximally possible. The values obtained in our research are gathered 
in Table Hm 

B. Systems comparison 

The main feature of a quantum circuit that gives us some information about the generated state is its size. In our 
work we put the main interest in the minimal number of quantum gates needed to obtain a circuit which manages to 
generate a maximally entangled state. The number of basic gates needed to generate maximum entanglement provide 
insight into the difficulty of this process. 

All obtained data are gathered in Tables [I] and |TT] The computation performed by us shows that the size of the 
circuit grows exponentially with the size of a system (see Figure [2]) . If we treat the gates used in our work as the 
elementary operations, we find out that the complexity of all quantum algorithms that harness the generation of 
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FIG. 2: The minimum number of gates (log scale) necessary for generating maximally (highly) entangled states. 



multipartite entanglement is exponential. It means that all hypothetical quantum algorithms which are supposed to 
bring essential computational speed-up by utilising multipartite entanglement may be practically inefficient. 

Additionally, considering both complete systems and spin chains, the obtained results show that the simulation 
of dynamics of a complete system using spin chain system brings an exponential growth of the number of quantum 
gates needed. That may suggest that quantum systems of informatics based on spin chains may be exponentially less 
efficient than the systems where the interactions between all pairs of qubits are possible. 



C. Optimal structures 



In our work we additionally analyse the structure of quantum systems for the presented circuits. While spin 
chain systems seem to be inefficient in the context of generating MME we find it interesting to find a structure that 
guarantees maximum efficiency without complete connections. The results shown in Figure [3] represent the structures 
of systems allowing the generation of maximum entanglement in an optimal time. These examples show that it is not 
necessary to allow the interaction between all qubits to obtain maximum efficiency. We can get an essential reduction 
of qubit connections without increasing the number of gates necessary to generate maximum entanglement. The 
evolution of the structure seems to be regular. 



IV. CONCLUDING REMARKS 



In this paper we have presented the method of applying the genetic algorithm to generate quantum circuits preparing 
maximally entangled states. The algorithm provided in this work can be applied to optimize the size of a circuit 
preparing maximally multipartite entangled states. We have provided an analysis of minimum number of quantum 
gates necessary for generating MMES. The size of a circuit is interesting both in the context of minimal complexity 
of a circuit generating the maximum entanglement and in obtaining maximally entangled states itself, because the 
complexity of the algebraic representation of a resulting state increases with the number of gates in a circuit. The 
circuit we have obtained for a 6-qubit system that contains 12 quantum gates and generates a maximally entangled 
state with 16 non-zero coefficients is an example. The best result known so far included a 13-gate circuit generating 
a state with 32 coefficients [217] . 

Moreover, the described method allows taking into account the topology of connections between the particles in the 
system. Developed method enables one to ensure that the resultant circuit will be possible to implement in a system 
with arbitrarily defined connections between particles. This gives us the possibility of analysis the process of generation 
of the multipatrtite entanglement. We exploit this possibility and provide the comparison of efficiency of generating 
maximally multipartite entangled states in the spin chain system and in the system of complete connections. We 




FIG. 3: Connection graphs of the optimal (a) 5-qubit (b) 6-qubit and (c) 7-qubit quantum circuit preparing 

maximally entangled state. 



have shown that in the case of generation the multipartite entanglement a spin chain is essentially less efficient than 
the system of complete connections. Moreover, we have presented some topologies with relatively sparse connections 
between qubits that enable the generation of maximum entanglement with the same efficiency as the systems of 
complete connections. 
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